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Electron momentum density and Compton profiles in Lithium along (100), (110), and (111) di- 
rections are calculated using Full-Potential Linear Augmented Plane Wave basis within generalized 



X 



gradient approximation. The profiles have been corrected for correlations with Lam-Platzman for- 
mulation using self-consistent charge density. The first and second derivatives of Compton profiles 
are studied to investigate the Fermi surface breaks. Decent agreement is observed between recent 
experimental and our calculated values. Our values for the derivatives are found to be in better 
agreement with experiments than earlier theoretical results. Two-photon momentum density and 



' O ■ one- and two-dimensional angular correlation of positron annihilation radiation are also calculated 

O , within the same formalism and including the electron-positron enhancement factor. 

> ' I. INTRODUCTION 

ON 

£ — . Compton scattering and positron annihilation techniques are well established tools for studying the momentum 
distribution of electrons in solids§~@. In the Compton scattering technique, the intensity distribution of energy broad- 

Os . 

ened Compton scattered radiation (called the Compton profile) is studied while in positron annihilation technique, 
measurement of angular correlation between two photons emitted during the annihilation of a thermalized positron 
with electrons in the solid, is performed. The Compton profile (CP) and angular correlation of positron annihilation 

o : 

radiation (ACPAR) curves contain the finger prints of the Fermi surface (FS) breaks in the momentum distribution in 



the first and higher Brillouin zones. It is well known that although the magnitude of the discontinuity in momentum 



distribution itself changes due to electron-electron, electron-ion and electron-positron correlations, the position of 
the discontinuity remains unchanged@~0. Thus, these techniques together are useful to extract information about FS 
geometry and to identify electron correlation effects. Positron annihilation techniques are more sensitive to the outer, 
weakly bound conduction electrons and are also capable of performing measurements in both one-dimensional and 
two-dimensional geometries, with much superior momentum resolution. Theoretically, however, it is more straight for- 
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ward to calculate CP's as ACPAR studies necessitate to account for electron-positron many-body correlation effects. 
These effects are incorporated in the form of momentum, energy or density dependent enhancement factor 

Early measurements of Compton profile were suffering from the limited momentum resolution in the Compton 
scattering experiment (~ 0.4a. it.). The advent of high intensity, high energy and well-polarized synchrotron sources 
and the spectrometers with high resolution (~ 0.12a.u.) have resulted in revival of interest in this area@~0. On the 
theoretical side, high performance computing facilities have made it possible to perform calculations on a fine p-mesh 
with a better convergence criteria for the total energy and charge density self-consistency. 

In the present communication, we report full potential linearized augmented plane wave (FP-LAPW) calculations 
of the electron momentum distribution in Li. In lithium, because of its smallness, the electron-ion interactions 
are strong and the electrons do not behave like a text-book example of homogeneous electron gas as in sodium. 
The Fermi surface of Li shows a small but definite departure from free electron sphere. The electron momentum 
density is highly anisotropic, the high momentum components (HMCs) are small but important. The momentum 
distribution of Li has been investigated in the past, both theoretically and experimentally, by several worker 
with the available state-of-the-art procedures. Recently, Sakurai et aZ.E^ have performed high resolution Compton 
scattering experiments for Li to measure Compton profiles and Fermi radii. The measured CP's are compared with 
the theoretical ones calculated using KKR method. Subsequently, Schiilke and coworkersEl measured 11— directional 
Compton profiles and employed them to reconstruct the three dimensional EMD in Li metal. The comparison of 
the experimental CP's with their theoretical counterparts have always shown some discrepancies. This is partly due 
to the various approximations involved in computing the profile and partly due to the experimental errors. The 
methods like FP-LAPW can provide accurate results for the one-electron wavefunction and hence should be employed 
for computation of theoretical CP's. Such accurate calculations of the CP's of Li have been reported by KuboEMl 
employing the GW approximation using FP-LAPW wavefunctions within local density approximation (LDA) as the 
zeroth approximation. The calculated CP's are in good agreement with the experimental values but the derivatives of 
the CP's match the experimental ones poorly0. The derivatives of CP's are important as they provide information 
about the Fermi surface. 

In the present paper, we use full-potential linearized augmented plane wave method for the computation of CP 
which have been corrected to include the correlation effects not accounted for within LDA, alongwith one and two 
dimensional angular correlation of positron annihilation radiation (ID- and 2D- ACPAR) which include corrections 
due to electron-positron correlation effects. Such studies using the same formalism for calculation of CP and ACPAR 
can provide complementary information about the Fermi surface and the electron momentum distribution. Our aim 
in this communication is to extract the FS geometry from the two techniques together and to identify the part of 
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electron correlations left out in the theory to describe the EMD. The plan of the paper is as follows: in section II, 
we present the momentum space formulation of the LAPW wavefunction and computational details, while section III 
deals with the results. 



II. METHODOLOGY AND COMPUTATIONAL DETAILS 

The one electron wavefunction for an electron in the state labeled by wave vector k and band index j is expanded 
in the Linearized Augmented Plane Wave (LAPW) basis </>g +( g(r) as 

G 

where c^ + g are the expansion coefficients and G denotes reciprocal lattice vector. The LAPW function is a plane 
wave outside the MT sphere while inside it is a linear combination of ui(r), the solution of the radial Schrodinger 
equation and its energy derivative ui(r). This allows greater flexibility inside the spheres (than the APW method) 
and hence permits computation of an accurate solutionS. 

The momentum space LAPW wavefunction is obtained by a Dirac- Fourier transformation of eq. (^) and the electron 
momentum density is computed from the momentum space wavefunction of occupied states. The Compton profile is 
then obtained by performing a double integral. 

The theoretical Compton profile thus obtained is overestimated at low momenta and is underestimated at higher 
momentum values than its experimental counterpart. This discrepancy is often attributed to the correlation effects 
that are ignored in the independent particle mode@. Lam and PlatzmanEl have shown that these effects can be 
incorporated into the EMD by augmenting to the independent electron model momentum density a correction term. 
The recent formulation of Cardwell and CooperEl based on the Lam-Platzman work which takes care of the non-unity 
occupation below the Fermi momentum kf and non-zero occupation beyond kf has been employed in the present 
work. 

For a single positron in a defect free crystal, the positron density will be distributed over the entire crystal. Therefore, 
the positron wavefunction, tp + (r), is considered to be delocalized and is represented by a plane wave basis. Further, 
as the positron essentially thcrmalizcs before annihilation, it is assumed to be in the state k + = 0. The enhancement 
in the TPMD due to the electron-positron short-range correlation is calculated according to the prescription given by 
Puska et aL0 which is a function of electron and positron densities. Thus, TPMD is evaluated using the following 
prescription : 



9^ 



OCC occ „ 

(P) = E I ! 2 = E I J ^■ r ^ + (^l(rWg(n e (rhn+(r)) d 3 r | 2 . (2) 
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where n e (f) and n + (f) are the electron and positron densities respectively and the enhancement factor g(n e (r), n + (r)) 
in the limit n + (r) — > is 

9a{r s ) = 1 + 1.23r s + 0.98890rf /2 - 1.4820r 2 + 0.3956r^ 2 + r 3 s /6 (3) 

,1/3 

where r. 
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A self-consistent band-structure calculation was performed using the LAPW method as implemented in WIEN97 
packageil. The calculation employs full potential which implies that the nonspherical part of the potential inside the 
muffin-tin sphere and its deviation from the constant potential in the interstitial region are taken into consideration. 
The simplified generalized gradient approximation (GGA) due to Perdew et aL0 was used for the exchange-correlation 
part of the Kohn-Sham potential. In the band structure calculation the lattice constant of Li in bcc structure is 
taken to be 6. 61375a. ufl The self-consistency cycles were carried out to an energy tolerance of 10~ 6 Ryd and 
charge convergence of 10~ 5 electrons. The various band-structure parameters agree very well with earlier accurate 
calculations!^!!]. 

For the calculation of Compton profile, in l/48th of the Brillouin zone (BZ) we used 40425 fc-points to evaluate 
the momentum space wavefunction which when translated by the reciprocal space vectors yield 40425 x 531 p-points. 
Linear tetrahedron mcthod0 was used for the calculation of Compton profiles over a momentum mesh of O.OOla.u.. 
The correlation correction was carried out as described by Cardwell et aZ@. In carrying out these calculations we 
used the self-consistent density inside the MT sphere whereas in the interstitial region we assumed the density to be 
flat with little structure. No r s cut-off value was used as suggested by Cardwell et a/El. 

For the calculation of two photon momentum density, the wavefunction for the positron was obtained by solving 
the secular determinant once using the self-consistent Coulombic potential from earlier calculation for electronic band 
structure but with an opposite sign. We have used 819 fc-points to evaluate Fl(p). 



III. RESULTS AND DISCUSSION 

The electron momentum density for lithium metal along < 100 >, < 110 > and < 111 > directions in the 
momentum space is plotted in Fig. 1. Each panel shows EMD along one direction, calculated within LDA and 
GGA employing von Barth-Hcdin (VBH)0 and Perdew-Burke-Ernzerhof (PBE)0 exchange-correlation potentials 
respectively alongwith GGA-EMD corrected using Lam-Platzman (LP)@ correction term in that direction employing 
the model momentum density proposed by Cardwell et a/.El. The results are almost identical for VBH and PBE 
exchange-correlation potentials indicating that the non-local corrections as described within GGA do not seem to 
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affect the momentum density in Li although the electronic structure is slightly different in LDA and GGA. This gives 
rise to almost identical Compton profiles (CPs) for LDA and GGA formulations as predicted by Lam and PlatzmanS. 
These observations for Li support the fact that although the ionic potential is strong in Li, the conduction electron 
density behaves more like a homogeneous electron gas. Similar calculations for transition metals do show significant 
differences in LDA and GGA resultsil. 

The strong ionic periodic potential, however, does couple the states near Brillouin zone (BZ) boundaries and the 
conduction electron wavefunctions contain strong high momentum components and the EMD does show anisotropic 
behavior as is evident from Fig. 1. Since Li has only one electron in the conduction band, in the one-electron picture, 
the EMD is zero beyond kf in the first BZ and between its images in higher zones. Lam-Platzman@ correction 
partly takes care of the correlation effects on the wavefunctions in one-electron picture. The occupation number of 
an interacting homogeneous electron gas is estimated to be smaller than that of a non-interacting free electron gas by 
4% for Li@. Our LP corrected EMD displays the effect of states below kf being pushed above kf in Li and these can 
be compared with results in Fig. 5 of ref.0 which uses mean occupation numbers derived from electron gas data for 
correlation effects and couples them with orthogonalized plane wave band structure method. EMD dies off in higher 
zones along < 100 > and < 111 > directions but shows a strong umklapp component along < 110 > direction. These 
higher momentum components (HMCs) play an important role in determining the shapes of Compton profiles and 
angular correlation curves. 

Compton profiles J^(q) corrected for the correlations are presented in Fig. 2 together with the experimental and 
KKR resultsEl with k along < 100 >, < 110 > and < 111 > directions. In order to facilitate the comparison 
with experimental CP's, the theoretical CP's were convoluted with a Gaussian of FWHM of 0.12 a.u. (experimental 
resolution in ref.El). Our Compton profiles, when compared with experimental values, support the widely known 
behavior, namely, overestimation at low momentum and underestimation at higher momentum values. Correlation 
corrections lower the CP values within the main Fermi surface but they are still higher than the experimental results. 
It is to be noted that KKR values are consistently higher than our values near q = and show a more pronounced 
cusp like behavior near kf for all the three directions. The cusp seen in CP reflects the discontinuity in EMD at kf 
in the first BZ and their images in higher zones. The results indicate that the discontinuity is smaller in FP-LAPW 
calculations than in KKR. The disagreement between present results and KKR results could be attributed to the fact 
that present work is a full potential calculation whereas KKR uses muffin-tin shape approximation. Secondly, the LP 
correction has been calculated using the prescription of Cardwell et aZ.El employing self-consistent charge density in 
the present work whereas KKR uses interpolated results from homogeneous electron gas data. Our results agree with 
those of KKR beyond q = 0.8 a.u. while the experimental values are consistently higher than the theory in this region. 
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This reflects that the non-zero occupation beyond kf for the inhomogeneous electron gas is only partially accounted 
for by the theory. 

We have also compared the directional anisotropics in the Compton profiles with the experimental work of Sakurai 
et a/Jll and the overall agreement is found to be good. The directional anisotropies are important while comparing 
theory with experiment as the systematic errors in the experimental and theoretical results are canceled out. The 
prominent structures near q = are well reproduced at the correct momentum values, however, they are overestimated 
by theory (Fig. 3). This is again due to the correlation correction functional as discussed by Bauer et a/.El. In the 
present calculation, we have included the correlation correction which is isotropic and therefore does not affect the 
results of anisotropy. Attempts to include anisotropic corrections are in progress and will be published elsewhere. 

The two-dimensional ACPAR surfaces were obtained by integrating p 2l (p) along < 100 >, < 110 > and < 111 > 
directions respectively; while the one dimensional curves were obtained using the linear tetrahedron method. The ID 
ACPAR curves along < 100 >, < 110 > and < 111 > are shown in Fig. 4 and are in qualitative agreement with 
earlier published wor The ID ACPAR curves were convoluted with a Gaussian with FWHM of 0.022 a.u.. The 

2D ACPAR data convoluted with 0.5x0.23 mrad 2 FWHmEEI is presented in Fig. 5 and the overall shape matches 
with experimentS. The higher momentum components in TPMD are seen quite clearly in the 2D ACPAR plot. The 
smaller bump comes from the TPMD in the second zone in the < 110 > direction, whereas the contributions from 
higher zones are hardly visible. The larger bump is a result of projected contributions from < 011 > and < 101 > for 
the same p value. The HMC intensity is seen to be a sharp function of momentum values. The momentum density 
decreases with increasing momentum values which is a direct reflection of s-like wavefunction. 

We have extracted relevant FS data from both electron and electron-positron momentum distributions. Although, 
in principle, both the Compton scattering and the ACPAR probe the electron momentum distribution and provide 
complementary information about FS, the latter however provides the best possible measurements for FS breaks 
in EMD since the electron-positron correlations enhance the momentum density at kf. The prominent breaks in 
momentum distribution are seen around 0.6 a.u. with slightly different values along different directions as shown in 
Fig. 1 for EMD. Similar structure is seen in TPMD also. The images of the FS breaks, seen in higher zones due 
to periodicity, are reflected in the CP and ACPAR data. The structures seen are identical to those seen earlier and 
discussed by Sakurai et a/Jll for CP and Stewart et a/.El for ACPAR curves. However, we point out that the LP 
correction shifts the FS breaks as seen in the derivatives of CP. This is a limitation of the model which employs 
isotropic momentum density to calculate the LP correction. 

Figure 6 displays the first derivatives of the directional Compton profiles. Although the breaks in the first derivative 
curve without convolution bring out the structures rather well indicating the distortions of the free electron sphere, 
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however, in Fig. 6, we show the derivatives of the convoluted Compton profiles to facilitate comparison with exper- 
imental results. The derivatives of ID ACPAR which bring out the FS breaks sharply are presented in Fig. 7. The 
values of the FS radii along the principal symmetry directions < 100 >, < 110 > and < 111 > as estimated from the 
positions of peaks in the second derivatives of CP differ from the actual Fermi radii as described by Sakurai et aZ.El. 
The first derivatives of directional ACPAR curves give the correct Fermi radii since the high momentum components 
in TPMD are smaller. Our values of Fermi radii alongwith the experimental valueJll are presented in Table I. The 
maximum Fermi surface asphericity, [fcno — kioo]/k® where k® is the free electron radius, in our calculation turns out 
to be 5.6% against the experimental value of 4.6%, 5% and 4.7% respectively obtained from CpEI, ID-ACPArEI and 
2D- ACPAR dataS. 

In this paper, we have presented the Compton profiles computed using FP-LAPW method within GGA and have 
corrected them for correlations on the lines of Cardwell and CooperEl using self-consistent charge density. The deriva- 
tives of the calculated CP's are in good agreement with their experimental counterparts. Usually, the discrepancy 
between theory and experiment is ascribed to the limitation of local density approximation. However, we have seen 
that gradient corrections to exchange-correlation potential as described in GGA do not affect the electron momen- 
tum density significantly. In principle, Lam-Platzman correction describes the non-local effects on the momentum 
density correctly but the practical implementation is able to account for it only partly, namely the isotropic elec- 
tron correlations are described satisfactorily but not the anisotropic ones. Although Bauer et aZ.El have rejected the 
idea of momentum density functional theory, we feel a description for the exchange-correlation energy functional in 
momentum space will allow a better quantitative representation of EMD theoretically. 

One and two dimensional ACPAR curves are also computed using the FP-LAPW method. We have shown that 
the different shapes of HMCs are well reproduced by our calculations. Inclusion of density dependent enhancement 
factor is found to reduce the HMCs at large p values. To our knowledge, this is the first theoretical report of ID- and 
2D-ACPAR for Li incorporating density dependent enhancement effects. 
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Table I 



Directions Fermi surface radii 

FP-LAPW Experimental 

fcmo 0.578 0.577 ± 0.004 

fcno 0.611 0.604 ± 0.004 

fcin 0.585 0.586 ± 0.004 
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FIG. 1. Electron momentum densities along < 100 >, < 110 > and < 111 > directions. 
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FIG. 2. Compton profiles along < 100 >, < 110 > and < 111 > directions. 
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FIG. 3. Anisotropies of the Compton profiles between various directions. 
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FIG. 4. One-dimensional ACPAR curves along < 100 >, < 110 > and < 111 > directions 
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FIG. 5. Two-dimensional ACPAR plot of Li in the (110) plane. 
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FIG. 6. Derivative of the Compton profiles along < 100 >, < 110 > and < 111 > directions. 
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FIG. 7. Derivative of one-dimensional ACPAR curves along < 100 >, < 110 > and < 111 > directions 



18 



